Printed on February 1, 2008 
Electronic and Structural Properties of Carbon Nano-Horns 
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We use parametrized linear combination of atomic orbitals calculations to determine the stabil- 
ity, optimum geometry and electronic properties of nanometer-sized capped graphitic cones, called 
"nano- horns" . Different nano-horn morphologies are considered, which differ in the relative location 
of the five terminating pentagons. Simulated scanning tunneling microscopy images of the various 
structures at different bias voltages reflect a net electron transfer towards the pentagon vertex sites. 
We find that the local density of states at the tip, observable by scanning tunneling spectroscopy, 
can be used to discriminate between different tip structures. Our molecular dynamics simulations 
indicate that disintegration of nano-horns at high temperatures starts in the highest-strain region 
near the tip. 
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Since their first discovery!]], carbon nanotubes have 
drawn the attention of both scientists and engineers 
due to the large, mimber of interesting new phenom- 
ena they exhibitaH, and due to .their potential use in 
nanoscakj devices: quantum wiresQ, nonlinear electronic 
elemental, transistors!!, molecular memory devicesB, and 
electron field emitterst3dl3. Even though nanotubes have 
not yet found commercially viable applications, projec- 
tions indicate that this should occur in the very near 
future, with the advent of molecular electronics and 
further miniaturization of micro-electromechanical de- 
vices (MEMS). Among the most unique features of nan- 
otubes are their electronic properties. I^-h&s been pre- 
dicted that single-wall carbon nanotubeai-lEj can be ei- 
ther metallic or semiconducting, depending on their di- 
ameter and chiraiity£Jj1i3. Recently, the correlation be- 
tween the chirality and conducting behavior of nanotubes 
has been confirmed by high-resolution scanning tunneling 
microscopy (STM) studiesNM 

Even though these studies haAje, demonstrated that 
atomic resolution can be achievedljEil, the precise de- 
termination of the atomic configuration, characterized 
by the chiral vector, diameter, distortion, and position 
of atomic defects, is still a very difficult task to achieve 
in nanotubes. Much of the difficulty arises from the fact 
that the electronic states at the Fermi level are only indi- 
rectly related to the atomic positions. Theoretical model- 
ing of STM images has been found crucial to correctly in- 
terpret experimental data for graphite 13 , and has been 
recently applied to carbon nanotubesc-JO As an alter- 
native technique, scanning tunneling spectroscopy com- 
bined with modeling has been used to investigate the 
effect of the terminating cap on the electronic structure 
of nanotubescl 

Among the more unusual systems that have been syn- 
thesized in the past. ff?w years are cone-shaped graphitic 
carbon structuresEZrcfl. Whereas similar structures have 
been obsepjed previously near the end of multi-wall 
nanotubesES, it is only recently that an unusually high 



production rate of up to 10 g/h has been achieved 
for single-walled cone-shaped structures, called "nano- 
horns" , using the CO2 laser ablation technique at room 
temperature in absence of a metal catalystETJ. These con- 
ical nano-horns have the unique opening angle of «20°. 

We consider a microscopic understanding of the elec- 
tronic and structural properties of nano-horns a cru- 
cial prerequisite for understanding the role of terminat- 
ing caps in the physical behavior of contacts between 
nanotube-based nano-devices. So far, neither nano-horns 
nor other cone-shaped structures have been investigated 
theoretically. In the following, we study the structural 
stability of the various tip morphologies, and the inter- 
relationship between the atomic arrangement and the 
electronic structure at the terminating cap, as well as 
the disintegration behavior of nano-horns at high tem- 
peratures. 

Cones can be formed by cutting a wedge from planar 
graphite and connecting the exposed edges in a seam- 
less manner. The opening angle of the wedge, called 
the disclination angle, is n(jr/3), with 0<n<6. This 
disclination angle is related to the opening angle of the 
cone by 6 — 2sin~ (1 — 71/6). Two-dimensional pla- 
nar structures (e.g. a graphene sheet) are associated 
with n = 0, and one-dimensional cylindrical structures, 
such as the nanotubes, are described by n — 6. All 
other possible graphitic cone structures with < n < 6 
have been obsecwed in a sample generated by pyroiysis 
of hydrocarbonsE3. According to Euler's rule, the termi- 
nating cap of a cone with the disclination angle n(n/3) 
contains n pentagon(s) that substitute for the hexagonal 
rings of planar graphite. 

The observed cone opening angle of «20°, correspond- 
ing to a 5ir/3 disclination, implies that all nano-horns 
contain exactly five pentagons near the tip. We classify 
the structure of nano-horns by distinguishing the rela- 
tive positions of the carbon pentagons at the apex which 
determine the morphology of the terminating cap. Our 
study will focus on the influence of the relative position 
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FIG. 1. Optimized carbon nano-horn structures with a 
total disclination angle of 5(n/3), containing five isolated pen- 
tagons at the terminating cap. Structures (a) - (c) contain 
all pentagons at the conical "shoulder", whereas structures 
(d) - (f) contain a pentagon at the apex. The pentagons are 
highlighted by a darker color. 

of these five pentagons on the properties of nano-horns. 

The cap morphologies investigated in this study are 
presented in Fig. |l|. Nano-horns with all five pentagons at 
the "shoulder" of the cone, yielding a blunt tip, are shown 
in Figs. |l](a)-(c). Nano-horns with a pentagon at the 
apex of the tip, surrounded by the other four pentagons 
at the shoulder, are shown in Figs. |l](d)-(f). Note that 
the cone angle of each nano-horn is «20°, even though 
the size of the terminating cap varies with the relative 
position of the pentagons. 

To determine the structural and electronic properties 
of carbon nano-horns, we used the parametrized linear 
combination of atomic orbitals (LCAO) technique with 
parameters determined by ab initio calculations for sim- 
pler structures^. This approach has been found useful 
to describe minute electronic structure and total energy 
differences for systems with too large unit cells to handle 
accurately by ab initio techniques. Some of the problems 
tackled successfully by this technique are the electronic 
structurC[-and superconducting properties of the doped 
Ceo solidES, the opening of pseudo^gaps near the Fermi 
level in a (10,10) nanotubes ropeEaMand a (5,5)@(10,10) 
double- wall nanotubeE3, as., well as fractional quantum 
conductance in nanotubesEJ. This technique, combined 
with the recursion technique to achieve an O(N) scaling, 
can determine very efficiently the forces on individual 
atomsED, and had previously been used with success to 
describe the disintegration dynamics of fullcrenesa, the 
growth of multi=wall nanotubestJ and the dynamics of a 
"bucky-shuttle"El. 

To investigate the structural stability and electronic 
properties of carbon nano-horns, we first optimized the 
structures with various cap morphologies, shown in 
Fig. |l| For the sake of an easier interpretation of our 
results, we distinguish the A^ cap w40 — 50 atoms at the 
terminating cap from those within the cone-shaped man- 
tle, that is terminated by N c a gc atoms at the other end. 



TABLE I. Structural data and stability results for 
carbon nano-horn structures (a)-(f), presented in Fig. |l]. 
A?tot = Ntip + iVcdgc is the total number of atoms, which 
are subdivided into tip and edge atoms. < i5 co h,tot > is the 
average binding energy, taken over the entire structure, and 
< ^coh.tip > the corresponding value excluding the edge re- 
gion. < i? C oh,edge > is the binding energy of the edge atoms, 
and < -Bcoh.pcnt > is the average over the pentagon sites in 
each system. 



Quantity 




(a) 


(b) 


(c) 


(d) 


(e) 


(f) 


iVtot 




205 


272 


296 


290 


308 


217 


N tip 




172 


233 


257 


251 


270 


180 


N c dgc 




33 


39 


39 


39 


38 


37 


< ^coh.tot 


> (eV) 


-7.28 


-7.29 


-7.30 


-7.30 


-7.31 


-7.28 


< -E co h,tip 


> (eV) 


-7.36 


-7.36 


-7.37 


-7.36 


-7.37 


-7.36 


< -Ecoh.cdgc > (eV) 


-6.88 


-6.88 


-6.88 


-6.88 


-6.87 


-6.89 


< Scoh.pcnt > (eV) 


-7.28 


-7.28 


-7.28 


-7.28 


-7.28 


-7.28 



We associate the tip region of a hypothetically infinite 
nano-horn with all the sites excluding the edge. Struc- 
tural details and the results of our stability calculations 
are presented in Table |. These results indicate that 
atoms in nano-horns are only wO.l eV less stable than 
in graphite. The relative differences in < i? C oh,tot > re- 
flect the strain energy changes induced by the different 
pentagon arrangements. To minimize the effect of under- 
coordinated atoms at the edge on the relative stabilities, 
we excluded the edge atoms from the average when cal- 
culating < i? C oh,tip >• Since our results for < i? C oh,tip > 
and < -Bcoh.tot > follow the same trends, we believe that 
the effect of edge atoms on the physical properties can 
be neglected for structures containing hundreds of atoms. 
Even though the energy differences may appear minute 
on a per-atom basis, they translate into few electron-volts 
when related to the entire structure. Our results suggest 
that the under-coordinated edge atoms are all less sta- 
ble than the cone mantle atoms by «0.5 eV. Also atoms 
in pentagons are less stable than those in hexagons by 
wO.l eV, resulting in an energy penalty of ~0.5 eV to 
create a pentagon if the strain energy induced by bend- 
ing the lattice could be ignored. 

When comparing the stabilities of the tip regions, de- 
scribed by < -Ecoh.tip >j we found no large difference be- 
tween blunt tips that have all the pentagons distributed 
along the cylinder mantle and pointed tips containing a 
pentagon at the apex. We found the structure shown 
in Fig. |l](c) to be more stable than the other blunt 
structures with no pentagon at the apex. Similarly, the 
structure shown in Fig. [l](e) is most stable among the 
pointed tips containing a pentagon at the apex. Equi- 
librium carbon-carbon bond lengths in the cap region 
are dec = 1-43 — 1.44 A at the pentagonal sites and 
dec — 1-39 A at the hexagonal sites, as compared to 
dec — 1-41 — 1.42 A in the mantle. This implies that the 
"single bonds" found in pentagons should be weaker than 
the "double bonds" connecting hexagonal sites, thus con- 
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FIG. 2. Simulated STM images in the tip region for 
(a) the nano-horn shown in Fig. ^(c) and (b) the nano-horn 
shown in Fig. |l|(d). These results for the occupied elec- 
tronic states near the Fermi level, corresponding to the bias 
voltage of Vt = 0.2 V, are suggestive of a net electron 
transfer from the hexagonal to the pentagonal sites. The 
charge density contour displayed corresponds to the value of 
p — I.35xl0 -3 electrons/A 3 . Dark lines depict the atomic 
bonds to guide the eye. 

firming our results in Table [I] and the analogous behavior 
in the Ceo molecule. 

Since pentagon sites are defects in an all-hexagon 
structure, they may carry a net charge. To character- 
ize the nature of the defect states associated with these 
sites, we calculated the electronic structure at the tip 
of the nano-horns. The charge density associated with 
states near Ep, corresponding to the local density of 
states at that particular position and energy, is propor- 
tional to the current observed in STM experiments. To 
compute the local charge density associated with a given 
eigenstate, we projected this state onto a local atomic 
basis. The projection coefficients were used in conjunc- 
tion with real-space atomic wave functions from density 
functional calculationsE3 to determine the charge density 
corresponding to a particular level or the total charge 
density. To mimic a large structure, we convoluted the 
discrete level spectrum by a Gaussian with a full-width 
at half-maximum of 0.3 eV. Using this convoluted spec- 
trum, we also determined the charge density associated 
with particular energy intervals corresponding to STM 
data for a given bias voltage. 

In Fig. ||, we present such simulated STM images for 
the nano-horns represented in Figs. |l](c) and 0(d). We 
show the charge density associated with occupied states 
within a narrow energy interval of 0.2 eV below the Fermi 
level!3 as three-dimensional charge density contours, for 
the density value of p — 1.35 x 10~ 3 electrons/A 3 . Very 
similar results to those presented in Fig. ^ were obtained 
at a higher bias voltage of 0.4 eV. As seen in Fig. ||, ppir 
interactions dominate the spectrum near Ep. These im- 
ages also show a net excess charge on the pentagonal sites 
as compared to the hexagonal sites. This extra negative 
charge at the apex should make pointed nano-horn struc- 
tures with a pentagon at the apex better candidates for 
field emittersE^rt 3 ] than structures with no pentagon at 
the apex and a relatively blunt tip. 

It has been shown previously that theoretical model- 
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FIG. 3. Local electronic densities of states, normalized 
per atom and shown in consistent arbitrary units, for the ter- 
minating cap of the tip (solid line) and for the pentagon sites 
only (dashed line) of nano-horn structures shown in Fig. 
The similarity between the results for the entire cap and the 
pentagon sites necessitated a vertical offset of one set of data 
for an easy distinction, thus indicating that variations in the 
arrangement of pentagons affect the densities of states not 
only at the pentagonal sites, but within the entire cap region. 
The most marked differences are noted near the Fermi level, 
at E^0 eV. 

ing of STM images is essential for the correct interpre- 
tation of experimental data. Atomically resolved STM 
images, however, are very harci-to obtain especially near 
the terminating caps of tubesc3 and cones, due to the 
large surface curvature that can not be probed efficiently 
using current cone-shaped STM tips. A better way to 
identify the tip structure may consist of scanning tunnel- 
ing spectroscopy (STS) measurements in the vicinity of 
the tip. This approach is based on the fact that in STS 
experiments, the normalized conductance (V/ I)(dl /dV) 
is proportional to the local density of states which, in 
turn, is structure sensitive. We have calculated the local 
density of states at the terminating cap of the tips for 
the different nano-horn structures shown in Fig. [y. Our 
results are shown in Fig. |^, convoluted using a Gaussian 
with a full-width at half-maximum of 0.3 eV. 

To investigate the effect of pentagonal sites on the elec- 
tronic structure at the tip, we first calculated the local 
density of states only at the 25 atoms contained in the 
five terminating pentagons. The corresponding densities 
of states, shown by the dashed curve in Fig. |^, are found 
to vary significantly from structure to structure near the 
Fermi levelcl Thus, a comparison between the densities 
of states at E^iEp should offer a new way to discriminate 
between the various tip morphologies. For an easy com- 
parison with experiments, we also calculated the local 
density of states in the entire terminating cap, including 
all five pentagons and consisting of iV«40 — 50 atoms, 
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depending on the structure. The corresponding density 
of states, given by the solid line in Fig. ||, is vertically dis- 
placed for easier comparison. Our results show that the 
densities of states, both normalized per atom, are very 
similar. Thus, we conclude that the pentagonal sites de- 
termine all essential features of the electronic structure 
near the Fermi level at the tip. ._. 

Next, we have studied the heat resilienceEj as well as 
the decay mechanism of nano-horns at extremely high 
temperatures using molecular dynamics simulations!!! 
In our canonical molecular dynamics simulations, we 
keep the structure at a_ constant temperature using a 
Nose-Hoover thermostats, and use a fifth-order Rungc- 
Kutta interpolation scheme to integrate the equations 
of motion, with a time step of At = 5 x 1CP 16 s. We 
found the system to remain structurally intact within 
the temperature range from T = 2, 000 — 4, 000 K. Then, 
we heated up the system gradually from T = 4, 000 K 
to 5, 000 K within 4,000 time steps, corresponding to a 
time interval of 2 ps. Our molecular simulations show 
that nano-horn structures are extremely heat resilient 
up to T<4, 500 K. At higher temperatures, we find these 
structures to disintegrate preferentially in the vicinity of 
the pentagon sites. A simultaneous disintegration of the 
nano-horn structures at the exposed edge, which also oc- 
curs in our simulations, is ignored as an artifact of finite- 
size systems. The preferential disintegration in the higher 
strain region near the pentagon sites, associated with a 
large local curvature, is one reason for the observation 
that nano-horn tips are opeaed easily at high tempera- 
tures, in presence of oxygenEl. 

In summary, we used parametrized linear combina- 
tion of atomic orbitals calculations to determine the 
stability, optimum geometry and electronic properties 
of nanometer-sized capped graphitic cones, called nano- 
horns. We considered different nano-horn morphologies 
that differ in the relative location of the five terminat- 
ing pentagons. We found a net electron transfer to the 
pentagonal sites of the cap. This extra charge is seen in 
simulated scanning tunneling microscopy images of the 
various structures at different bias voltages. We found 
that the local density of states at the tip, observable by 
scanning tunneling spectroscopy, can be used to discrim- 
inate between different tip structures. Our molecular dy- 
namics simulations indicate that disintegration of nano- 
horns at high temperatures starts in the highest-strain 
region near the tip. 

We thank S. Iijima for fruitful discussions on carbon 
nano-horns. We acknowledge financial support by the Of- 
fice of Naval Research and DARPA under Grant Number 
N00014-99-1-0252. 
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